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ABSTRACT 

Mazeika's  method  for  forecasting  mixed-layer  (thermocline)  depth 
of  the  upper  ocean  layers  is  discussed  along  with  a  newer  version  of 
this  method  developed  by  James.   Using  Mazeika's  method  primarily,  a 
verification  for  the  Northeast  Pacific  Ocean  was  completed  with  data 
from  Ocean  Weather  Stations  PAPA  (50N,  145W)  and  NOVEMBER  (30N,  140W) 
and  a  point  named  MIDPOINT  (40N,  140W) .   The  results  indicate  Mazeika's 
method  is  successful  at  Station  PAPA  more  than  seventy-five  percent  of 
the  time  during  the  heating  season  followed  by  a  rapid  decline  as  the 
cooling  season  begins.   The  method  should  be  useful  in  the  entire 
Central  Subarctic  Domain  as  described  by  John  P.  Tully.   The  method  fails 
at  NOVEMBER  and  MIDPOINT  producing  less  than  thirty  percent  success  in 
prediction.   James'  version  did  not  improve  the  results  obtained  at 
Station  NOVEMBER.   This  failure  appears  to  be  due  to  the  controlling 
parameters  for  processes  in  the  Subtropic  or  Transitional  oceanographic 
regions  (which  include  NOVEMBER  and  MIDPOINT);  these  differ  from  para- 
meters controlling  oceanic  processes  in  the  Pacific  Subarctic  region 
(Station  PAPA) ,  which  resemble  those  involved  in  the  Atlantic  region 
for  which  Mazeika's  method  was  developed. 

Climatology  data  which  can  be  used  to  obtain  surface  and  400-foot 
level  temperature  are  also  tested.   The  results  indicate  these  data  are 
very  useful  and  accurate  in  determining  the  stability  index  required  of 
Mazeika's  method. 
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1.   Background. 

The  United  States  Navy  is  expending  much  effort  in  the  direction 
of  anti-submarine  warfare.   The  proper  use  of  anti-submarine  detection 
devices  hinges  on  the  capability  of  correctly  predicting  the  thermo- 
cline  depth  in  the  various  oceans.   Any  reliable  method  which  can  be 
easily  and  properly  used  by  the  men  and  equipment  of  the  anti-submarine 
forces  would  add  to  the  overall  defense  capability  of  our  country.   To 
this  end,  the  testing  of  one  particular  method  of  thermocline  depth 
prediction  was  undertaken. 


2.   Introduction. 

Mazeika's  method  utilizes  a  subjective  but  realistic  approach  to 
the  thermocline  depth  prediction  problem.   The  thermocline  depth  pre- 
dicted by  this  procedure  is  the  mean  depth  of  the  interface  between 
the  mixed  layer  and  the  thermocline.   The  method  requires  the.  use  of 
easily  obtainable  inputs  such  as  forecast  wind  speed,  fetch,  wind 
duration  and  the  synoptic  weather  picture.   Subjective  aspects  of  this 
method  come  to  light  in  certain  choices  which  are  made  by  the  user  of 
the  prediction  technique,  as  will  be  explained  later. 

Various  types  of  vertical  mixing  are  included  in  Mazeika's  method. 
They  are  mechanical  mixing,  which  is  due  to  the  combined  effect  of  wind 
waves  and  associated  wind  currents;  instability  mixing  due  to  advection 
and  density  changes;  and  mixing  due  to  vertical  convergence -divergence 
caused  by  changes  in  the  overlying  atmospheric  pressure  field,  or  by 
convergence -divergence  in  the  surface  current  field.  Mechanical  mix- 
ing is  introduced  into  the  method  by  use  of  a  sea  state  parameter  which 
is  based  on  the  wind  wave  characteristics  of  the  area  under  investiga- 
tion.  The  use  of  a  stability  index  and  corrections  for  the  salinity 
gradient  were  chosen  to  satisfy  the  requirements  of  instability  mix- 
ing.  Finally,  convergence -divergence  is  accounted  for  in  the  construc- 
tion of  graphs  used  to  determine  the  mixed-layer  depth. 

Mazeika  combined  the  use  of  these  mixing  parameters  with  the  study 
of  BT  data,  weather  charts  and  simultaneous  weather  and  wave  observa- 
tions from  the  Atlantic  Ocean  Weather  Station  (OWS)  CHARLIE  (52N,  35W) 
to  develop  his  prediction  procedure  empirically.  After  development  was 
completed,  data  from  Atlantic  OWS  BRAVO  (56. 5N,  51W) ,  ECHO  (35N,  48W) , 
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and  DELTA  (44N,  41W)  were  used  to  test  the  method.   No  data  from  the 
Pacific  ocean  weather  stations  were  used  for  the  empirical  develop- 
ment of  the  method  or  for  its  testing. 

Mazeika  has  stated  that  his  prediction  method  can  generally  be 
used  anywhere  in  the  ocean  without  previous  information  on  the  mixed 
layer  and  thermocline;  the  only  requirement  is  knowledge  of  the  surface 
wind  waves  and  the  stability  index.   The  testing  of  Mazeika' s  method  in 
the  Pacific  Ocean  would  serve  to  substantiate  this  claim  by  Mazeika  or 
possibly  bring  to  light  weaknesses  in  the  method  which  may  be  present. 
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3.   Mazeika's  Method. 

Mazeika  (1960)  considered  the  mixed-layer  thickness,  h,  to  be  a 
function  of  wave  amplitude  A,  wavelength  ?\t    and  a  mixing  parameter  k, 
which  is  valid  at  the  interface  between  the  mixed  layer  and  the  thermo- 
cline.   The  mixing  parameter  cannot  be  expected  to  remain  constant. 
Since  it  applies  to  the  bottom  of  the  mixed  layer,  it  must  be  dependent 
on  mixed-layer  thickness  and  on  the  wave  parameters.   In  addition,  its 
value  depends  on  the  stability  in  the  thermocline. 

The  relation  describing  orbital  motion  in  trochoid  waves 

includes  most  of  the  parameters,  if  h  is  considered  to  be  a  fixed  value 
under  given  surface  conditions  and  if  stability  in  the  thermocline  is 
constant.   The  above  function,  k(A,h,  ?0 ,  was  adopted  by  Mazeika  in 
view  of  the  generally  accepted  theory  that  orbital  motion  of  particles 
due  to  waves  decreases  exponentially  with  depth  and  is  proportional  to 
the  wave  parameters. 

Solving  (1)  for  h, 


fc-&M#)-  <2> 


Equation  (2)  can  be  used  to  compute  h,  if  the  wave  parameters  and 
the  corresponding  value  of  k  for  existing  stability  in  the  thermocline 
are  known. 

In  the  above  equations  the  amplitude  (A)  is  one-half  of  the  sig- 
nificant wave  height  (H.  ._)  ;  wave  length  (/*\)  was  computed  by  the  re- 

o 

lation  ^=  3.41T   ,  where  T    is  the  period  of  maximum  energy  of  the 
max        max  w 

spectrum  for  a  fully-developed  sea. 
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On  the  assumption  that  vertical  temperature  difference  is  propor- 
tional to  the  mean  temperature  gradient  in  the  upper  part  of  the  thermo- 
cline,  Mazeika  adopted  the  difference  in  the  entire  thermocline  as  the 
index  of  stability  (a  salinity  correction  could  be  added,  see  Table  I). 
Since  most  BT  data  usually  reach  a   depth  of  400  feet,  Mazeika  selected 
the  difference  between  the  surface  and  the  400  foot  level  as  the  stabil- 
ity index 

25  -  t0  -  t400. 

A  deeper  layer  (e.g.,  600  feet)  could  be  taken  if  data  were  available. 

A  certain  initial  stability  exists  with  a  given  AT  in  the  absence 
of  a  mixed  layer.   If  the  mean  stability  index  remains  constant,  stabil- 
ity increases  from  the  initial  value  with  increasing  mixed-layer  thick- 
ness.  Increase  of  stability  is  slow  when  the  mixed  layer  begins  to 
form  and  gradually  becomes  greater  as  the  mixed-layer  thickness  and  AT 
increase. 

To  connect  mixing  parameter  values  with  stability  and  wave  para- 
meters, another  parameter  (sea  state  parameter)  was  introduced  by  Mazeika 


T\  "  Hi/3T, 


max 


fl  is  used  as  a  measure  of  sea  state  for  wind  waves,   k  values  deter- 
mined from  equation (1) with  known  wave  parameters  and  mean  mixed-layer 
thickness  (h)  values  were  plotted  against  T\  values  computed  with 
corresponding  wave  parameters  for  eight  groups  of  AT.   Curves  which  fit 
the  distribution  of  points  for  each  AT  group  were  computed  and  resulted 
in  the  solid-line  curves  of  figures  1-8.   h  values  which  were  related 
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to  the  k(T\)  curves  were  then  computed  and  plotted  as  the  dashed-line 
curves  on  figures  1-8.   The  central  h  curve  applies  to  normal  mixing, 
the  upper  h  curve  applies  to  convergent  situations  with  horizontal 
convergence  occurring  in  the  upper  ocean  layers  adjacent  to  the  point 
of  interest,  while  the  lower  h  curve  applies  to  situations  with  hori- 
zontal divergence.  A  later  discussion  covers  procedures  used  in  deter- 
mining which  h  curves  apply  for  a  particular  location.   The  h  curve  is 
used  to  determine  the  mixed-layer  thickness  in  feet  from  the  h  scale  on 
the  left  margin  of  figures  1-8. 

Richard  W.  James  (1966)  of  the  Naval  Oceanographic  Office  has  de- 
veloped a  new  version  of  thermocline  depth  prediction  based  entirely 
on  Mazeika's  method.   This  author  had  completed  most  of  the  Pacific 
verification  data  before  James'  version  was  published  so  that  only  lim- 
ited testing  of  it  was  introduced  in  this  report. 

James  believes  Mazeika's  model  to  be  the  most  realistic  approach 
to  wind  mixing  to  date,  but  points  out  that  no  account  is  taken  of  the 
demonstrated  time  lag  between  onset  of  wind  and  the  mixing  subsequently 
generated,  and  also  that  it  is  overly  sensitive  to  low  wave  conditions. 
James  has  corrected  these  deficiencies  in  new  figures  which  are  described 
below,  and  has  redefined  the  stability  index  to  what  he  believes  is  a 
more  realistic  parameter. 

James  retained  the  sea  state  parameter, TV,  as  defined  by  Mazeika, 
but  he  replaced  its  determination  method  by  developing  a  new  graph 
(figure  9)  which  is  based  on  wave  forecasting  procedures  of  Pierson, 
Neumann,  and  James  (1955).   This  can  be  used  in  all  cases  of  fully- or 
non-fully-developed  sea. 
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To  account  for  the  time  lag  between  occurrence  of  wind  and  subse- 
quent mixing  James  developed  Table  II.   This  table  gives  approximate 
time  periods  for  a  given  1\  value  to  stabilize,  since  a  certain  wind  re- 
quires a  given  time  to  develop  a  maximum  1\ value.   If  this  time  require- 
ment is  not  met  then  a  lower  1\  value  must  be  considered. 

James  also  believed  that  Mazeika's  stability  index,  AT,  is  a  fair- 
ly conservative  factor  and  that  it  may  not  be  representative  of  the  re- 
sistance to  mixing  in  many  cases.   Thus,  he  redefined  the  stability 
index  (AT..)  as  the  temperature  gradient,  in  °F/100  feet,  at  the  most 
shallow  thermocline.   The  thermal  gradient  does  not  have  to  extend 
through  100  feet  since  only  the  rate  of  change  is  desired. 

Salinity  corrections  for  this  method  are  obtained  from  figure  10. 
The  correction  is  given  in  °F/100  feet  and  is  added  to  the  previously 
determined  stability  index  (AT  ).   The  amount  of  correction  is  propor- 
tional to  the  salinity  gradient,  water  temperature,  and  the  mixing  fac- 
tors. 

James  also  developed  new  mixed-layer  depth  forecasting  curves  from 
the  formula: 


LD=^(i-tr^AT^ 


where  K  and  K  are  constants.   These  new  curves  are  depicted  in  figure 
11.   Convergence  and  divergence  was  accounted  for  in  this  figure  by 
using  the  next  AT.  curve  above  the  normal  AT-  curve  for  convergence  , 
and  the  AT.  curve  below  the  proper  value  for  divergence. 

The  effects  of  swell  on  the  mixing  of  water  was  also  considered 
by  James,  who  states  that  a  much  reduced  value  of  T\  should  be  used  for 
swell.   The  reduction  of  the  mixing  force  with  swell  is  due  to  the 
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narrowness  of  the  period  band,  lack  of  breaking  crests,  and  reduced 
angular  spreading.   He  also  states  that  swell  caused  from  decaying  wind 
waves  in  a  fetch  still  possesses  some  of  the  characteristics  of  sea  and 
should  represent  about  50  percent  of  f\  while  swell  arriving  from  a 
distant  storm  would  represent  approximately  25  percent  of  T\.      This  re- 
duction is  strictly  a  subjective  estimate. 
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4.   Mazeika's  Prediction  Procedure. 

To  use  the  method  properly,  the  stability  index  must  be  determined 
first  and  corrected  for  salinity  gradient  if  possible.  This  determina- 
tion can  be  accomplished  in  various  ways  based  on: 

(1)  averaging  temperature  differences  between  the  surface  and 
the  400-foot  level  for  four  to  six  bathythermograph  (BT)  observations 
taken  within  the  time  interval  of  24  hours,  preferably;  the  time  inter- 
val could  be  longer,  however,  up  to  a  maximum  of  ten  days. 

(2)  climatology  data  records  (monthly  or  seasonal)  which  give 
sea  surface  temperature  and  the  temperature  at  400  feet.  (The  accuracy 
of  these  data  will  be  discussed  later.) 

(3)  synoptic  sea  surface  chart  and  climatological  data  for 
400  feet. 

After  AT  has  been  determined,  the  salinity  gradient  correction 
can  be  applied  if  such  data  are  available.   The  salinity  correction  is 
found  by  entering  Table  I  with  the  mean  temperature  in  the  mixed  layer 
and  the  salinity  difference  (  /oo)  in  the  mixed  layer  and  coming  out 
with  a  value  AT'.   This  value  is  then  used  to  determine  the  correction 
C  =  -ST'-l*  for  negative  salinity  gradients  or  C  =2ST,-1°  for  positive 
salinity  gradients,  since  k(f\)  curves  actually  correspond  to  1°F  high- 
er stability  indexes.   This  C  correction  is  then  applied  to  the  origi- 
nal ZJf  value  to  determine  the  proper  graph  for  the  mixed  layer  predic- 
tion. 

The  next  step  in  the  prediction  procedure  is  to  determine  the  sea 
state  parameter  1\.   This  can  be  done  by  the  use  of  Table  III  which  lists 
values  of  wave  parameters  and  minimum  fetch  and  duration  for  various 
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windspeeds  for  a  fully-developed  sea.   For  application  in  a  non-fully- 
developed  sea,  H. /0  and  T    are  taken  from  graphs  in  H.O.  603  and  is 

l/J      max 

then  computed  with  the  equation  T\  =  H.  ._T 

The  forecast  windspeed  can  be  a  critical  factor  in  determining  an 
accurate  X\  •   The  error  in  mixed-layer  depth  prediction  caused  by  erro- 
neous windspeed  is  greater  with  a  small  stability  index  than  with  a 
larger  one.   To  illustrate,  an  error  of  1  knot  in  windspeed  could  pro- 
duce a  difference  of  16  feet  in  the  predicted  mixed-layer  depth  using 
a  stability  index  of  2°  while  this  same  error  would  produce  a  differ- 
ence of  6  feet  using  a  stability  index  of  13°.   The  accurate  forecast 
of  windspeed  is  imperative  to  an  accurate  mixed- layer  prediction. 

Finally,  a  decision  must  be  made  as  to  whether  an  area  is  under- 
going horizontal  divergence  or  convergence  which  may  influence  the  ef- 
fectiveness of  the  mixing  processes.   This  is  a  difficult  problem  which 
is  approached  subjectively  through  a  study  of  the  forecast  area.   Some 
areas  may  be  under  the  influence  of  permanent  convergent  or  divergent 
systems  which  may  preclude  conventional  mixing  processes.   To  aid  in 
this  determination  Mazeika  has  illustrated  in  figures  12  and  13  what 
he  believes  to  be  areas  of  strong  convergence  and  divergence.   Figures 
12  and  13  were  used  as  guide  lines  by  this  author  in  the  verification 
procedure  and  they  seemed  to  work  quite  well.   Laevastu  and  Hubert  (1965) 
note  that  convergence  and  divergence  usually  depend  on  a  subjective 
estimation  and  they  recommend  the  alternative  of  calculation  of  diver- 
gence and  convergence  from  surface  current  fields. 

After  the  required  parameters  have  been  determined  and  the  appro- 
priate graph  selected  the  thermocline  depth  can  be  determined.   The  fl 
value  is  found  on  the  far  left  margin  of  the  graph  (figures  1-8) . 
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From  the  horizontal  intersection  of  this  value  with  the  proper  k( f\) 
solid-line  curve  (normal,  convergent  or  divergent),  proceed  vertically 
to  intersect  the  corresponding  dashed  curve,  h(k).   Then  proceed  hori- 
zontally to  the  h  scale  on  the  left  margin  to  determine  the  thermocline 
depth. 
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5.   Verification  Procedure. 

The  requirements  for  the  use  of  Mazeika's  method  are  explicit. 
To  meet  these  requirements  a  search  for  data  was  initiated  and  led  im- 
mediately to  the  elimination  of  all  data  which  included  only  one  BT  per 
day.   These  sources  were  mainly  commercial  and  military  ship  reports. 
Since  data  from  Ocean  Weather  Stations  PAPA  (50N,  145W)  and  NOVEMBER 
(30N,  140W)  were  readily  available  at  the  Naval  Postgraduate  School,  a 
decision  was  made  to  limit  the  verification  to  the  Northeast  Pacific 
Ocean.   Also,  the  difference  in  geographical  location  of  these  two 
stations  in  the  Northeast  Pacific  Ocean  area  would  serve  to  test  the 
procedure  in  perhaps  different  ocean  regions  within  this  broad  area. 

Data  for  Station  PAPA  were  taken  from  booklets  prepared  by  the 
Canadian  Oceanographic  Data  Center.   This  source  is  considered  excel- 
lent as  all  parameters  were  included  in  it  except  the  synoptic  weather 
picture.   Ocean  Station  NOVEMBER  BT  data  were  obtained  from  the  Navy 
Fleet  Numerical  Weather  Facility  (FNWF)  at  Monterey9  California.   These 
data  were  in  the  form  of  daily  reports  which  normally  included  from  two 
to  ten  BT  observations. 

A  point  between  these  two  ocean  weather  stations,  which  will  be 
referred  to  as  MIDPOINT  (40N,  140W) ,  was  also  included.  Information 
for  MIDPOINT,  furnished  by  FNWF,  was  available  for  only  one  month. 

The  synoptic  weather  charts  which  were  needed  to  estimate  conver- 
gent, divergent  or  normal  synoptic  situations  were  obtained  from  records 
kept  in  the  meteorology  storeroom  at  the  Naval  Postgraduate  School  and 
also  from  the  Synoptic  Weather  Maps  (Daily  Series)  available  through 
the  U.  S.  Government  Printing  Office.  Most  of  the  charts  used  were 
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facsimile  charts  sent  from  the  National  Meteorological  Center  at 
Washington,  D.  C.   The  particular  synoptic  chart  used  included  surface 
pressure  and  1000-500  MB  thickness  analysis  for  the  Northern  Hemisphere. 

All  data  collected  for  the  verification  were  selected  from  months 
during  which  the  seasonal  thermocline  is  a  prominent  feature.   May 
through  August  was  considered  the  best  period  for  this  to  occur  at  all 
three  locations.   Station  PAPA  data  included  continuous  observations 
from  May  1964  through  October  1964.   September  and  October  were  included  to 
determine  the  method's  effectiveness  after  the  seasonal  thermocline 
period  has  ended.   A  limited  amount  of  observations  for  June  1966  were 
included  also.   The  data  periods  for  Station  NOVEMBER  included  the  fol- 
lowing: July  1957,  July  through  September  1958,  May  and  June  1964,  June 
1965  and  April  through  June  1966.   Data  for  MIDPOINT  were  available  for 
May  1965. 

As  Mazeika  points  out  in  his  report,  a  large  amount  of  continuous 
data  (which  must  include  all  necessary  parameters)  from  the  same  observa- 
tion point  is  not  easy  to  obtain.   The  data  gathered  for  this  verifica- 
tion were  found  in  many  different  places  and  brought  together  so  that 
all  requirements  were  met.   Numerous  other  continuous  data  sources  had 
to  be  discarded  as  one  or  two  parameters  could  not  be  located  for  the 
time  period  involved. 

After  all  available  data  had  been  collected,  a  verification  form 
was  developed  to  include  all  the  required  information  necessary  to  make 
proper  hindcasts.   This  form  is  presented  in  figure  14.  Most  of  the 
column  headings  are  self-explanatory,  but  explanations  for  a  few  are 
required. 

When  available,  windspeed  was  obtained  from  weather  observations 
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taken  concurrently  with  the  BT  observations;  otherwise,  it  was  esti- 
mated from  the  daily  synoptic  weather  charts  and  entered  on  the  form. 
Much  effort  was  expended  in  making  correct  windspeed  estimates  as  large 
errors  could  develop.  The  magnitude  of  this  type  of  error  was  brought 
out  earlier  in  this  report.  After  windspeed  was  determined,  Table  III 
was  entered  to  obtain  the  value  T\ . 

Next,  observed  mixed- layer  depth  was  entered  as  the  average  obtain- 
ed from  all  BT  observations  for  that  day.   Since  Mazeika's  method  pre- 
dicts the  mean  mixed- layer  depth  it  must  be  compared  with  an  average 
value  (since  individual  BT  observations  may  differ  greatly).   Tempera- 
tures at  the  surface  and  400  feet  were  averaged  and  then  compared  to 
obtain  an  accurate  stability  index. 

After  all  parameters  were  obtained,  the  proper  AT  graph  was  used 
in  selecting  the  mixed-layer  thickness  (thermocline  depth).  This  value 
was  then  compared  with  the  observed  mixed- layer  depth.   If  the  observed 
value  were  larger  than  the  predicted  one,  the  difference  between  the 
two  was  entered  as  a  negative  number  in  the  Difference  column.   If  the 
predicted  value  were  the  larger  of  the  two  values,  a  positive  number 
was  entered  in  the  Difference  column.   This  was  done  to  determine  if 
Mazeika's  method  was  over-  or  under-predicting  the  mixed-layer  thick- 
ness in  the  Pacific  Ocean. 

When  the  proper  criteria  have  been  met,  a  check  mark  is  placed  in 
the  column  labelled  Verification.   In  this  study  a  tolerance  of  +20 
feet  from  the  observed  value  is  considered  to  have  verified.   The 
tolerance  limits  were  chosen  after  considering  the  factors  described  in 
the  following  paragraph. 

First,  Mazeika  has  shown  in  figure  15  that  the  average  amplitude 
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of  internal  waves  is  17.85  feet.   This  is  interpreted  to  indicate  that 
internal  waves  alone  could  provide  a  variation  of  thermocline  depth  of 
this  size  about  its  mean  value.   Second,  Tabata  and  Giovando  (1962) 
have  shown  this  same  magnitude  in  figure  16.   They  also  concluded  that 
this  variation  was  due  to  internal  waves.   Finally,  any  method  which 
can  correctly  predict  within  20  feet  of  the  established  thermocline 
depth  is  certainly  a  valuable  tool  for  the  operating  ASW  forces. 
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6.   Verification  Results. 

To  insure  that  Mazeika's  method  was  being  applied  correctly  so 
that  the  results  obtained  would  be  useful,  a  trial  run  on  Atlantic  OWS 
CHARLIE  was  initiated.  A  period  of  one  month  was  tested  with  the  results 
presented  in  Table  IV.  The  method  predicted  within  20  feet  of  the  ob- 
served mixed-layer  depth  63.6%  of  the  time.   This  was  considered  accept- 
able, since  the  hindcasts  were  made  with  no  previous  experience  on  new 
data,  although  in  the  ocean  region  for  which  the  method  was  developed. 
After  this  trial  run,  it  was  believed  the  method  was  being  applied  cor- 
rectly. 

Station  PAPA  had  a  total  of  106  hindcasts  performed  during  the  heat- 
ing period  May  through  August.   Of  these,  82  were  verified  by  the  ob- 
served layer  depth  to  show  77.4%  of  the  hindcasts  verified  within  20 
feet  of  the  observed  average  layer  depth  of  that  day  (Table  V).  In- 
cluding those  hindcasts  which  verified  within  30  feet  of  the  observed 
value  increased  the  verification  to  85%.   Since  the  mean  observed  thermo- 
cline  depth  was  approximately  100  feet  (Table  VI),  the  above  percentages 
can  be  directly  related  to  verification  data  presented  by  Mazeika  (1960). 
These  percentages  are  slightly  higher  than  those  which  Mazeika  observed 
on  data  from  Stations  BRAVO  and  DELTA  in  the  Atlantic.   However,  all 
results  obtained  from  Station  PAPA  compared  favorably  with  what  Mazeika 
had  presented. 

A  monthly  analysis  of  Station  PAPA  verification  results  presented 
in  Tables  VII,  VIII,  and  IX  shows  the  high  success  of  the  method  during 
the  heating  season.  Also  shown  is  a  rapid  decline  in  the  method's  ef- 
fectiveness as  the  cooling  season  begins.   Table  VIII  indicates  46%  suc- 
cess of  prediction  was  experienced  during  the  month  of  September  and 
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only  21%  during  October.   The  period  of  seasonal  or  "transient"  thermo- 
cline  has  ended  and  the  mixed  layer  begins  to  extend  to  its  permanent 
thermocline  depth;  this  prediction  method,  which  depends  on  wind  and 
mechanical  mixing,  ceases  to  be  effective  since  it  does  not  include  a 
way  of  accounting  for  convective  instability  in  the  ocean  mixed  layer. 

Station  PAPA  is  located  in  a  position  similar  in  some  respects  to 
the  position  of  Station  CHARLIE,  at  which  Mazeika  originally  developed 
his  method.   The  weather  is  quite  similar  in  that  many  pressure  systems, 
which  normally  consist  of  high  winds  and  associated  moderate  to  severe 
storms,  travel  eastward  through  this  area.   The  mixing  of  the  upper  lay- 
ers of  water  due  to  these  disturbances  is  probably  very  similar  at  both 
locations.  Also,  the  latitude  of  Stations  CHARLIE  and  PAPA  is  very 
nearly  the  same.   This  would  indicate  that  the  temperature  and  radia- 
tion conditions  are  similar.   Acknowledging  the  fact  that  water  masses 
do  vary  greatly,  the  above  similarities  would  lead  one  to  a  general  ac- 
ceptance that  Mazeika' s  method  should  be  expected  to  work  quite  well  at 
Station  PAPA. 

Returning  to  Table  V,  it  can  be  seen  that  Station  NOVEMBER  dis- 
plays a  completely  different  set  of  results  than  that  of  Station  PAPA. 
NOVEMBER  verified  only  52  out  of  191  attempts  for  a  verification  score 
of  27.2%.   Table  VI  displays  a  mean  negative  value  of  error  much  higher 
than  the  mean  positive  value.   This  indicates  that  the  method  was  pre- 
dicting a  mixed-layer  depth  much  too  shallow  for  the  NOVEMBER  area. 
Station  NOVEMBER  is  in  a  semi -convergent  mixing  zone,  so  that  the  area 
is  normally,  but  not  continually,  subjected  to  convergent  mixing  pro- 
cesses.  The  NOVEMBER  area  is  usually  under  the  influence  of  high  pres- 
sure systems  which  move  eastward.   The  isobar  spacing  in  these  pressure 
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systems  is  usually  quite  large  resulting  in  low  windspeeds.   The  low 
windspeed  produces  a  reduced  value  of  sea  state  parameter  which  ulti- 
mately leads  to  a  shallower  forecast  of  mixed-layer  thickness,  using 
Mazeika's  method. 

Thus  the  NOVEMBER  area  is  quite  different  from  the  development 
area  of  Mazeika's  method  in  many  respects.   The  temperature  and  radia- 
tion conditions  differ  as  well  as  the  windspeed  regimes  and  the  overall 
weather  patterns.   These  differences  are  believed  to  be  the  reasons  for 
the  failure  of  Mazeika's  method  in  this  area. 

With  the  failure  of  Mazeika's  method  clearly  evident,  it  was  thought 
a  test  using  James'  version  would  be  in  order  to  see  if  the  new  stability 
index  and  forecasting  curves  might  improi^e  the  results  at  OWS  NOVEMBER. 
As  seen  in  Table  IV  the  results  were  also  very  poor.   Only  one  month's 
data, previously  tested  using  Mazeika's  method,  were  used  for  this  trial 
run. 

Even  with  the  poor  results  displayed,  two  interesting  side  effects 
were  noted.   First  of  all,  obtaining  the  temperature  gradient  of  the 
shallowest  thermocline  did  not  prove  to  be  an  easy  task.   From  the  data 
given  in  various  BT  reports,  the  gradient  is  determined  as  only  a  sub- 
jective estimate  at  best.   There  were  some  times,  however,  when  fairly 
accurate  gradients  could  be  determined.   These  occurred  when  sufficient 
information  was  reported  in  the  thermocline  gradient  zone.   Secondly, 
the  forecast  curves  of  James  are  a  much  improved  version  compared  to 
those  of  Mazeika.   They  are  very  easy  to  interpolate,  especially  at 
low  sea  state  parameter  values. 

Since  both  Mazeika's  method  and  James"  version  did  not  prove  to 


be  successful  at  Station  NOVEMBER,  one  might  attempt  to  improve  the 
forecasting  there  through  a  new  set  of  curves  empirically  developed  as 
was  done  at  Station  CHARLIE.   This  could  also  lead  to  a  set  of  curves 
which  would  be  useful  at  lower  latitudes  in  general  where  similar  condi- 
tions prevail. 

In  the  hope  of  locating  a  transition  area  where  the  method  works 
somewhat  better  than  at  NOVEMBER,  a  point  approximately  half-way  be- 
tween NOVEMBER  and  PAPA  was  tested.  MIDPOINT  had  27  hindcasts  perform- 
ed with  8  verifications  for  a  29.6%  predictability  figure  (Table  V). 
This  figure  indicated  that  one  must  be  somewhere  above  40N  latitude  in 
order  for  the  method  to  work  more  effectively.  MIDPOINT  also  had  a 
large  negative  error  value  indicating  that  the  method  was  still  not 
forecasting  deep  enough  thermocline  depths. 

Tully  (1964)  describes  the  behavior  of  the  Northeast  Pacific  Ocean 
and  suggests  through  his  description  where  the  transition  between  the 
NOVEMBER  regime  and  PAPA  regime  may  occur.   In  his  paper,  Tully  des- 
cribed various  domains  or  regions  which  require  various  models  for 
determination  of  temperature  structure  and  behavior.  These  domains 
are  shown  in  figure  17.  As  can  be  seen  MIDPOINT  lies  in  the  lower 
part  of  the  Transitional  Domain  and  Station  PAPA  lies  in  the  center  of 
the  Central  Subarctic  Domain.   Each  domain  is  considered  to  have  in- 
dividually appropriate  meteorological  and  oceanographic  mechanisms 
working  within  it.   If  this  is  true  and  the  study  by  Tully  indicates 
this  to  be  so,  then  it  can  be  assumed  that  Mazeika's  method  should  work 
effectively  within  the  Central  Subarctic  Domain.   This  is  not  meant  to 
limit  the  use  of  Mazeika's  method  but  only  to  point  out  an  area  where 
high  predictability  reasonably  may  be  expected. 
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In  summary,  Mazeika's  method  does  produce  effective  results  in 
the  Station  PAPA  area  and  perhaps  these  results  may  be  extrapolated  to 
include  the  whole  Central  Subarctic  Domain.  At  about  45N  the  method's 
forecast  ability  may  begin  to  lessen.  At  MIDPOINT  and  Station  NOVEMBER 
the  method  is  less  than  30%  reliable,  indicating  the  need  for  another 
consideration  of  the  forecasting  problem  in  this  area;  James'  version 
of  the  forecast  method  did  not  improve  results  at  NOVEMBER.  Mazeika's 
method  was  easy  to  apply  in  most  respects  except  for  low  values  of  T\^ 
where  curves  are  difficult  to  interpolate  accurately.   Some  errors  in 
forecast  mixed-layer  thickness  may  have  occurred  due  to  this  problem. 

James'  version  uses  a  stability  index  which  may  be  difficult  to 
determine  in  some  cases.   The  thermocline  gradient  is  not  always  clear- 
ly reported  and  thus  must  be  estimated.  The  overall  procedure  is  easy 
to  use,  especially  with  regard  to  James'  forecasting  curves  which  offer 
many  improvements  over  the  earlier  Mazeika  version. 
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7.    Verification  of  Climatology  Results. 

For  the  determination  of  stability  index  Mazeika  suggested  as  a 
source  climatology  data  (monthly  or  seasonal)  which  gives  the  tempera- 
ture at  the  surface  and  400  feet.   Since  actual  BT  data  were  used  for 
the  determination  of  stability  index  in  this  report,  a  comparison  of 
the  two  was  considered  for  a  reliability  check  of  the  climatology  data. 

Only  certain  areas  have  been  studied  thoroughly  enough  to  pro- 
duce such  data.   One  of  these  areas  is  the  Gulf  of  Alaska  in  the  North- 
east Pacific  Ocean.   Margaret  K.  Robinson  (1957)  of  Scripps  Institution 
of  Oceanography  conducted  a  study  of  sea  temperature  for  this  area. 
Robinson  analyzed  16,103  individual  BT  observations  from  the  period 
1941  to  1952.   The  BT  data  were  subjected  to  extrapolation,  interpola- 
tion and  visual  smoothing  to  obtain  the  various  monthly  averages. 

Figures  18  and  19  show  one  version  of  the  temperature  averaging 
results.   These  charts  can  be  used  to  obtain  a  monthly  average  stabil- 
ity index.   However,  a  different  version  may  be  more  accurate  for  deter- 
mining temperature  from  a  day-to-day  or  week-to-week  basis.   This 
version  is  displayed  in  figure  20.   This  graphical  summary  was  used  to 
obtain  the  temperature  data  for  the  comparison. 

Temperatures  for  the  surface  and  400  feet  were  taken  from  figure 
20  for  the  months  of  May  through  October.   These  were  compared  against 
BT  observation  averages  which  resulted  from  this  verification  report. 
The  results  have  been  tabulated  and  are  presented  in  Table  X. 

It  must  be  pointed  out  that  figure  20  is  a  summary  of  data  for  49N, 
148W  which  is  slightly  southwest  of  Station  PAPA.   All  temperatures 
taken  from  figure  20  were  slightly  higher  than  those  which  might  be 
expected  at  PAPA,  but  a  comparison  should  still  prove  valuable. 
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Since  data  from  Station  PAPA  were  available  for  the  last  half  of 
May  only,  Table  X  includes  the  averages  for  this  period  only.   It  can 
be  seen  that  the  averages  compare  quite  well  and  in  most  cases  would 
result  in  using  the  same  forecasting  curve  for  Mazeika's  method.   These 
averages  ought  to  be  even  closer  if  both  were  done  for  the  same  observa- 
tion point. 

From  this  comparison  it  is  evident  that  climatology  data  can  be 
used  effectively  in  determining  accurate  stability  indexes  in  Mazeika's 
thermocline  depth  forecasts. 
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8.    Conclusions  and  Recommendations. 

Mazeika's  method  of  thermocline  depth  prediction  does  work  ef- 
fectively at  Station  PAPA.   The  method  predicted  within  20  feet  of  the 
observed  mixed-layer  depth  77.47.  of  the  time.   This  is  considered  to 
be  a  good  prediction  percentage. 

From  the  study  conducted  by  Tully,  it  can  be  assumed  the  same 
success  experienced  at  Station  PAPA  may  be  realized  for  the  entire 
Central  Subarctic  Domain.   This  is  due  to  the  consistent  meteorologi- 
cal and  oceanographic  properties  of  this  area. 

Station  NOVEMBER  offered  little  success  with  either  the  Mazeika 
or  James  version.   Both  predicted  within  20  feet  of  the  thermocline 
depth  less  than  307.  of  the  time. 

The  verification  attempt  at  MIDPOINT  resulted  in  less  than  307. 
prediction  success  also.   The  success  of  Mazeika's  method  at  PAPA  and 
the  failure  of  it  at  NOVEMBER  and  MIDPOINT  is  primarily  due  to  weather 
and  environmental  conditions.   The  temperature,  radiation  conditions, 
wind  regime  and  weather  similarities  with  regard  to  Station  PAPA  and 
the  development  area  of  Mazeika's  method  leads  one  to  believe  the  method 
should  work  under  these  similar  conditions.   At  Station  NOVEMBER  and 
MIDPOINT  the  conditions  are  not  the  same  thus  causing  the  method  to 
fail.   It  is  believed  that  forecasting  curves  could  be  developed  for 
this  area  which  would  make  the  method  useful  here. 

During  the  verification  process  it  was  found  that  Mazeika's  fore- 
casting curves  were  difficult  to  interpolate  at  low  sea  state  para- 
meter values.   This  difficulty  may  lead  to  errors  in  mixed-layer  depth 
prediction.   All  other  parameters  and  procedures  were  considered  very 
easy  to  use.   This  method  should  be  adaptable  to  fleet-wide  distribu- 
tion for  use  in  anti-submarine  operations. 
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James9  version  of  Mazeika's  method  was  also  easy  to  use,  with  the 
exception  of  the  stability  index.   James"  thermocline  gradient  is  more 
difficult  to  determine  than  the  temperature  difference  which  Mazeika 
uses.   The  forecasting  curves  developed  by  James  are  much  easier  to 
use  and  provide  for  a  more  accurate  mixed-layer  forecast. 

Climatology  data  which  may  be  used  to  determine  the  stability 
index  for  Mazeika's  procedure  were  also  tested.   The  results  indicated 
that  climatology  data  supply  accurate  information  for  this  determina- 
tion. 
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Figure  12     'RELATIONSHIP  OF  TYPICAL  ASSUMED  CONVERGENT 
AND  DIVERGENT  FIELDS  AS  APPLIED  TO  PRESSURE 
SYSTEMS  FOR  CHOICE  OF  tj   CURVES 
(C  =  CONVERGENCE,     D=  DIVERGENCE) 

(Mazeikas    196, 
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Figure  13      RELATIONSHIP  OF  TYPICAL  ASSUMED  CONVERGENT 
AND  DIVERGENT   FIELDS  AS  APPLIED  TO  PRESSURE 
SYSTEMS  FOR  CHOICE  OF  77  CURVES 
(C  =  CONVERGENCE,        D=  DIVERGENCE) 
(Mazeikas    1960) 
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